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The IIB matrix model proposes a mechanism for dynamically generating four dimensional space- 
time in string theory by spontaneous breaking of the ten dimensional rotational symmetry SO (10). 
Calculations using the Gaussian expansion method (GEM) lend support to this conjecture. We 
study a simple SO(4) invariant matrix model using Monte Carlo simulations and we confirm that 
its rotational symmetry breaks down, showing that lower dimensional configurations dominate 
the path integral. The model has a strong complex action problem and the calculations were made 
possible by the use of the factorization method on the density of states p„(x) of properly nor- 
malized eigenvalues X„ of the space-time moment of inertia tensor. We study scaling properties 
of the factorized terms of p n (x) and we find them in agreement with simple scaling arguments. 
These can be used in the finite size scaling extrapolation and in the study of the region of con- 
figuration space obscured by the large fluctuations of the phase. The computed values of X n are 
in reasonable agreement with GEM calculations and a numerical method for comparing the free 
energy of the corresponding ansatze is proposed and tested. 
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1. Introduction 

Matrix models have been studied intensively in the past few years in the context of non- 
perturbative formulations of string theory and in the study of gauge/gravity duality By dimen- 
sionally reducing D = 10 dimensional \J(N) supersymmetric (SUSY) Yang-Mills theories to zero 
dimensions, one obtains the IIB Matrix Model ^ (IKKT model) which has been proposed as a 
non-perturbative definition of IIB superstring theory. In this model space-time is represented by 
the distribution of eigenvalues of the bosonic matrices, a feature that raises the possibility of dy- 
namical compactification of the extra dimensions by Spontaneous Symmetry Breaking (SSB) of 
the SO(D) rotational symmetry of the model. Such a scenario is plausible as calculations using the 
Gaussian Expansion Method (GEM) indicate [Q]. 

Monte Carlo simulations of matrix models [Q, could play an important role in understanding 
string theories in a similar fashion that lattice QCD has contributed to the understanding of the 
non-perturbative regime of quantum field theories. Unfortunately such simulations are plagued 
by the complex action problem which arises when one simulates the system after integrating out 
the fermionic degrees of freedom. This problem is particularly important in the lattice studies of 
finite density QCD [||]. The factorization method has been proposed in [^] as a general method to 
reduce the complex action problem and eliminate the overlap problem, see also [f7|]. The basic idea 
is to control an appropriately chosen variable in order to sample regions of the configuration space 
which are hard to sample using reweighting and whose contribution is crucial in the computation of 
the physical observables. The study of the scaling properties of the related density of states allow 
for useful extrapolations to the physical results. 

We present preliminary results from calculations performed on a related zero-dimensional 
matrix model proposed in ^ which realizes the scenario of dynamical compactification of space- 
time dimensions [|L Hj. The model has a very strong complex action problem and bears strong 
similarities to the IIB matrix model, which makes it a useful playground for testing ideas to apply 
on the IIB matrix model and more generally on other interesting physical systems with a complex 
action problem. We are able to show that SSB occurs consistently with the predictions in [||, g]. The 
scaling properties of the density of states are studied in detail and are found to agree with simple 
scaling arguments. This is possible only by sampling heavily suppressed regions by using the 
factorization method and it is crucial in the extrapolations used in order to compute the expectation 
values of the SSB order parameters. 

2. The Model 

Consider the partition function ^ 

Z = J dAdydye-( Sb+s ^ where S b = ^tr(A M ) 2 , 5 f = -xj/ f a (F^^A^ . (2.1) 

An (jU = 1 , . . . , D, D even) are. N xN hermitian matrices, and and Ya w& A^-dimensional row 
and column vectors.The actions 5b and Sf have an S\J(N) symmetry. The spinor index a = 1, . . . ,p, 
where p represents the number of components of a D-dimensional Weyl spinor, p = lPl 2 ~ x , and 
the flavor index /=!,... , Aff, where Aff represents the number of flavors. The p x p matrices are 
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Figure 1: The VEV (A„)o (n — 1,2,3 and 4) in the phase quenched model Zq are plotted for r = 1 (left) and 
r = 2 (right) against The data for each r can be nicely fitted to straight lines meeting at the same point 
(1 + r/2) atN = °°, which demonstrates the absence of SSB. 
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Figure 2: The small-jc (left) and large-jc (right) behavior of — r (x) for r = I. The straight lines are fits to 
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the theoretical behavior (3.1). 
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Figure 4: (Left) The function -^logw„(x) is plotted together with the scaling function Q n {x) extracted 



from the asymptotic behaviors (3.2) for r = 1 and n = 2. (Right) The solution to (2.5) is obtained by finding 

1 (o) , d 

the intersections of —~fn (x) and — <t>„ (x) for r = I and n = 2. The position of the intersections is indicated 

N dx 

by the arrows with the symbols x s and x\ for the regions x < 1 and x > 1, respectively. 



SO(D) gamma matrices after the Weyl projection. Thus the actions (gJJ) have an SO(D) symmetry, 
where the bosonic variables transform as vectors and the fermionic variables transform as Weyl 
spinors. Integrating out the fermions, we obtain Z = / dAe~ Sb Zf[A], where Zf[A] = (det^)^ and 
& = r^A^ is a pN x /W matrix. The fermion determinant detiF for a single flavor is complex 
in general. Under parity transformation Ad — > —Ad, A,- — > A,- (i 7^ D), the fermion determinant 
transforms as det & — > (det This implies that det & is real for configurations with Ad = and 
that the phase of the determinant becomes stationary for configurations with Ad = Ad-\ =0. We 
take the large-A^ limit with r = Nf/N fixed, which corresponds to the Veneziano limit. Whether 
the SSB of SO(D) occurs in that limit is the issue we would like to address. For that purpose, 
we consider the "moment of inertia tensor" T^ v = ^tr(A^A v ) and its real positive eigenvalues X n 
(n = 1, . .. ,D) ordered as X\ > A2 > • • • > Xd- The vacuum expectation values (VEV) of these 
eigenvalues (A„) play the role of the order parameters. If they turn out to be unequal in the large-Af 
limit, it signals the SSB of SO(D). Consider the D = 4 case where we have p = 2 and = <7,-, 
i = 1,2,3 and T^ = 104. The "phase-quenched model" is defined by 



jdAe- So [A] , S [A] = 5b [A] — Nf log | det ®\A] \ 



It is easy to show that the absence of SSB implies 



1 + ^ for all n = 1,2,3,4, 



(2.2) 



(2-3) 



where the VEVs ( • }o are taken with respect to fl2.2|) , which is confirmed at infinitesimal r [§], and 
also at r = 1 and r = 2 numerically in this work. In the full model, GEM calculation up to 9-th 
order [^] indicate that the true vacuum is only SO(2) invariant and the (A„) are not all equal. 

In order to simulate (2.1) we rewrite it as Z = J dAe.~ So ^ e' r W. This system is very hard to 



simulate using simple reweighting due to the complex action and overlap problem which make 
the simulations of large systems exponentially hard. In this work we use the factorization method 
proposed in ref. [H] where one computes the density of states of a properly chosen observable 
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Table 1: The solutions (x s , x\) to eq. ( 2.4 ) that correspond to the (local) maxima of p„ (x) are shown. We also 
add the corresponding VEV (A„) obtained by the Gaussian expansion method in |^]: xsO(2) = (Ki) obtained 
using the SO(2) ansatz and xso(3) = (Ki) obtained using the SO(3) ansatz. Bold/italic numbers for same r 
and n are to be compared according to the text. 



by studying a set of systems where the observable is constrained to a given fixed value. The 
choice of the eigenvalues A„ is promising since restricting their values to be large or small favors 
configurations with relatively small fluctuations of the phase T. In this respect it is convenient to 
define \i=jj^ an d the density of states p„(x) = (^8(x — and p^\x) = (^8(x — /L„)^> . Then 

it is easy to show that p n (x) = ^ p„\x)w n (x), where C={e' r )o = (cosr)o- It follows that (A„) = 
Jq dxx p n (x) and the deviation of its value from one is a measure of the effect of the phase. The 
function w n (x) is defined by w n (x)=(e' T }, lrX = (cosr)„ a , where ( • } na denotes a VEV with respect 
to the partition function Z„, v = / dAe~ So 8{x — X n ). It turns out that w„(x) > 0, which simplifies 
our analysis significantly. Using the saddle point approximation, the problem of determining (Xn) 
can be reduced to that of minimizing the "free energy" J^ n (x) = — logp„(x) by solving the saddle 
point equation 

4~ log Pn (X) = /i 0) (X) + 4- log W n (jc) = , (2.4) 

dx dx 

where /i (x) = ^logp„ (x). It is important that the errors due to statistics and finite do not 
propagate exponentially to (X n ) as a direct computation would imply. 

The implementation of the above system is obtained by studying Z n y = / dAe^^ So+v ^\ 
where V(z) = ^y(z — and 7 and E, are real parameters. The parameter 7 controls the position 
and width of the peak of \ n and it is chosen large enough so that the results become independent 
of its value. In our simulations we used 7 in the range 10 3 -10 7 . Using the fact that p n y{x) = 
^5(x — Xi)^ p^ (x) exp I — V'^x(A 7! }o^ j, where ( • ) n y is a VEV with respect to Z n y, the 
position of the peak of the distribution function p„y(x) is given by the solution of 

/i 0) (x) - (A„)oV'(x(A„)o) = . (2.5) 

If we denote the solution by x p , we use the estimators x p = (A„)„y, w n (x p ) = (cosr)„y and 

/i 0) (x p ) = (W((Xn)ny) = r(An> {(K)n.y - 

3. Results 

First we study the phase quenched model Zq. We simulate the system for r = 1,2 and compute 
the eigenvalues (A„)o. We find that (X n )o{N) = 1 + r/2 + @{l/N) as can be seen from fig. [j] We 



5 



A study of the complex action problem 



Konstantinos N. Anagnostopoulos 



conclude that no SSB occurs in the phase quenched model and verify eq. (23). We calculate f„°^ (x) 



from eq. (^5|) by simulating Z n V - Using simple scaling arguments we find that its asymptotic 
behavior at x <C 1 and x 3> 1 is 

1 f^M~ / {h( 5 - n ) + r5 n\)\ + a n x<l , 

at 2 '" w "\-W^>o+(5+0; *>!■ 

Similar arguments lead to the respective asymptotic behavior of w„(x) 

1, / \ / \ I -c n x 5 -" (x< l,n = 2,3,4) , 
^lnw„(x)^(x) = _^_ (4 _ n) (X>>M=1)2)3) . M 



By varying the constants a n , c„ and d„, we fit our data to eqs. ( |3.1[ ) and (p3|). We verify the 
expected asymptotic behaviors and use the coefficients c„ and d n in order to extrapolate <!>„ (jc) to 
the region in x where we find the solution to the saddle point equation (2.5). In fig. ^| we show the 



scaling ( |Q| ) for r = 1 and in fig. || the scaling Q_2) for r = 1 and n = 2. The solution to eq. (2.5) 



is determined from the intersection of the curves ^fji {*) with — & n (x) for each n. The results for 



r = 1 and « = 2 are shown in fig. [| 

We use the notation x s and x\ for the solutions in the x < 1 and x > 1 regions respectively 
which correspond to the local maxima of p„(x). For n = 1 we obtain only x/ = (Ai) and for n = 4 
we obtain only x s = (A4). We tabulate the results in table [l] and we compare them with those 
obtained using GEM in [^]. We note that since the dominant configurations near x\ for n = 2 and x s 
for n = 3 are typically two dimensional, these are to be compared with the SO(2) ansatz. Similarly, 
since the dominant configurations near x/ for n = 3 and x s for n = 4 are typically three dimensional, 
these are to be compared with the SO(3) ansatz. 

We find that (X\) > 1 > (A4), a relation that is clearly going to survive the large-Af limit. 
Therefore we conclude that SO(4) SSB manifests in the model. In order to determine the group 
that SO (4) breaks to, we need to calculate the dominant peak as ,/V — > °°. We consider the quantity 
A n = -jj { log p n (x/ ) - log p n (x s ) I = 3>„ (xi ) - ® n (x s ) + S n , where E„ = ^ dx | ^/i 0) (x) } . If A„ > 
the peak at x\ dominates, otherwise x s . We find A2 ~ 0.34 for r = 1 and A2 ~ 0.25 for r = 2 and we 
conclude that SSB breaks at least down to SO (2). Unfortunately A3 turns out to be very close to 0, 
so we are unable to determine if SO (4) breaks to SO (2) as GEM predicts or to SO (3). 



4. Conclusions 



We have tested a scenario for dynamical compactification of space-time by simulating a toy 
matrix model related to the IIB matrix model of string theory. We have shown SSB of SO(4) 
rotational symmetry consistent with GEM analysis and small r calculations[||, §]. The phase 
quenched model has no SSB, confirming the expectation that the wild fluctuations of the phase 
of the fermionic partition function plays a crucial role in the mechanism of SSB. Large and small 
length scales are dynamically generated by these fluctuations which make the calculation of the 
dominant ones in the thermodynamic limit a challenging problem. Our results indicate how to 
proceed with the study of the IIB matrix model. Although the latter case is computationally more 
demanding, SUSY could make SSB easier to see. 
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Calculations were possible because of the use of the factorization method. By effectively 
sampling large and small x regions we are able to exploit the asymptotic behaviors of pi°^ (jc) and 
w n {x) in order to extrapolate the results to regions in x and system size which are inaccessible by 
direct simulations of the phase quenched model. It is possible that a remaining overlap problem 
makes the results of table [l] slightly differ from GEM results. By a generalization of the factor- 
ization method[ 10] this problem can be overcome and achieve also better quantitative agreement. 
Then Monte Carlo studies of many interesting systems hindered by the complex action problem 
are hopefully going to be made possible by using the factorization method. 
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